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Abstract 



We analyze 9.7 x 10 6 -B-B pairs recorded with the CLEO detector to deter- 
mine the production ratio of charged to neutral S-meson pairs produced at 
the T(4S) resonance. We measure the rates for B° -» J/ipK^ and B+ -> 
J/x^K^ + decays and use the world-average S-meson lifetime ratio to extract 
the relative widths fy^ = r (T(4SHB+Er) = L04 ± o.07(stat) ± OM(syst). 

Too r(T(4S)— >B°B ) V ' V ' 

With the assumption that + /oo = 1, we obtain /oo = 0.49 ± 0.02(stat) ± 

O.Ol(syst) and f-\ = 0.51±0.02(stot)±0.01(syst). Lhis production ratio and 

its uncertainty apply to all exclusive I?-meson branching fractions measured 
at the T(4S) resonance. 
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Measurements of exclusive S-decay branching fractions from e + e~ collider operation at 
the T(4S) resonance assume equal production rates of charged and neutral 5-meson pairs [0. 
In the literature, the uncertainty in a specific branching fraction due to a lack of knowledge 
of the production ratio is often ignored. 

Any physics based upon comparisons of absolute decay rates of charged and neutral B 
mesons will profit from a more precise knowledge of the 5-production ratio, /+-//00 = 
r(T(4S') — > B + B~)/T(T(4S) — > B°B ). For example, a comparison of the branching 
fractions of two-body hadronic decays can be used to obtain information on the relative 
contributions from external and internal spectator decays 0. For all exclusive decay modes 
studied, the B + branching fraction was found to be larger than the corresponding B° branch- 
ing fraction, indicating constructive interference between the external and internal spectator 
amplitudes. This is in contrast to the destructive interference observed in hadronic charm 
decay. The magnitude of the constructively interfering fraction depends on the value of 
f+-/ foo- Another application of the / + _// o ratio arises in the use of ratios of charmless 
hadronic B-decay rates to set bounds on the angle 7, the phase of the CKM matrix 
element V u b |],|J. The uncertainty on / + _// o contributes to the systematic uncertainty of 
the 7 bound. 

A better measurement of / + _//oo would also allow a more meaningful comparison with 
theoretical predictions of the relative B + B~ and production rates at the T(4S) 

resonance. If there are no other important differences between the two T(4S) decays, 
such as B + — B° mass splitting or isospin-violating form factors in the decay amplitude, 
Coulomb corrections to B + B~ production near threshold are not negligible, giving rise to 
r(r(AS)-*B+B ) ^ |gn othgj. au thors [El argue that the S-meson substructure cannot 

be ignored and strongly reduces the Coulomb effect in the S-production ratio to 1.05—1.07, 
depending on the B masses and momenta. 

Existing measurements of the admixture ratio of charged to neutral B mesons produced at 
the T(4S) resonance have an uncertainty of ~15%. One measurement |7| used the branching- 
fraction ratio of B(B + -> J/ipK^ + ) to B(B° -> J/ipK^ ) f§ to yield ^ x ^ = 1.15 ± 
0.17±0.06, where the first uncertainty is statistical, the second is systematic, and tb denotes 
the B lifetime. Another measurement |§ used a ratio of B — > D*lv decays to extract 
h=L x ^+ = 1.14 ±0.14 ±0.13. 

700 T g 

In the present analysis, we study the decays B — > J/tpK^*\ which are isospin conserving 
transitions, since the J/ip daughter is an iso-singlet and the B and mesons are both 
iso-doublets. The decays B + — > J/^pK^ + and B° — > J/i/jK^ must therefore have equal 
partial widths and we can extract R = j^- x = ^jf^u~/^^(*)o^ ; where M is the efficiency- 
corrected signal yield. Using the ratio of two similar decay rates to extract R, we exploit 
the cancellation of common experimental uncertainties. Throughout this Letter, reference 
to charge conjugate states is implicit. 

The data analyzed in this study were recorded at the Cornell Electron Storage Ring 
(CESR) with two configurations of the CLEO detector, CLEO II and CLEO II. V. The data 
consist of an integrated luminosity of 9.2 fb _1 of e + e~ annihilations recorded at the T(4S') 
resonance and of 4.6 fb -1 taken in the continuum, 60 MeV below the T(4S) energy. The 
results in this Letter are based upon 9.7 x 10 6 BB candidates and supersede those of 
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Ref. [0. 

The components of the CLEO detector most relevant to this analysis are the charged- 
particle tracking system, the 7800-crystal Csl electromagnetic calorimeter, and the muon 



chambers. The first third of the data were collected with the CLEO II detector ||10|| , which 
measured the momenta of charged particles in a tracking system consisting of an inner 6-layer 
straw-tube chamber, a 10-layer precision drift chamber, and a 51-layer main drift chamber, 
all operating inside a 1.5 T solenoidal magnet. The main drift chamber also provided a 
measurement of the specific ionization loss (dE/dx) used in particle identification. Two 
thirds of the data were taken with the CLEO II. V configuration, for which the innermost 



straw-tube chamber was replaced with a 3-layer silicon vertex detector ||11|| , and the argon- 
ethane gas of the main drift chamber was replaced with a helium-propane mixture. The 
muon identification system in both the CLEO II and CLEO II. V configurations consisted of 
proportional counters placed at various depths in the return yoke of the magnet. 

Since the backgrounds for B — > J/ipK^*' decays are very low, track and photon quality 
requirements have been designed to maximize signal yield. We reconstruct B — > J/^K^ 
candidates in the data samples taken at the T(4S') energy. 

Candidate J/ip mesons are reconstructed in their leptonic decay modes, requiring J/ip 
lepton daughter tracks to have momenta larger than 800 MeV/c. 

For J/ip reconstruction in the muon channel, one of the muon candidates was required to 
penetrate the steel absorber to a depth greater than three nuclear interaction lengths. For 
the opposite sign daughter candidate, no muon detection requirement was imposed. 

Electron candidates were identified based on the ratio of the track momentum to the 
associated shower energy in the Csl calorimeter and specific ionization loss in the drift 
chamber. Bremsstrahlung produces a radiative tail in the e + e~ invariant mass distribution 
below the J/tj) pole. We recovered some of the resultant efficiency loss by detecting the 
radiated photon. We selected photon candidates (E y > 10 MeV) with the smallest angle 
to the e ± track, provided this angle did not exceed 5°. The J/ip — > e + e~ efficiency was 
increased by ~20%, without adding background. 

We reconstructed 15 900 ± 700 inclusive J/ip — > l + l~ candidates (Fig. [I]), about equally 
shared in the two dilepton reconstruction modes. The resolution in the J/tp invariant mass 
was ~13 MeV. 

We required the dimuon invariant mass to be within 50 MeV of the world-average J/ip 
mass 0, corresponding to a ~3.5 standard deviation (a) selection. For the dielectron in- 
variant mass we required —150 MeV < (m ee — mj/^) < 50 MeV to allow for the radiative 
tail. The J/ip energy resolution was improved by a factor ~4 after performing a kinematic 
fit of the dilepton invariant mass to the J/ip mass. We required J/ip candidates to have 
momenta below 2 GeV/c, which is near the kinematic limit for J/ip mesons originating from 
a B meson nearly at rest. 

The Kg — > 7r + 7r~ candidates were selected from pairs of tracks forming well-measured 
displaced vertices. The resolution in 7r + 7r~ invariant mass is approximately 2.5 MeV. Due to 
very low background in B —>■ J / ipK s candidates, we only require that neutral kaon candidates 
have a normalized mass within 10cr (because the Kg mass distribution has non-negligible 
non-Gaussian tails) and a normalized flight distance greater than zero. 

Charged kaon and pion candidates are required to have a measured dE/dx within 3cr 
of the energy loss expected for the given particle type. Neutral pions are reconstructed 
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FIG. 1. Invariant mass spectrum of the (a) J ftp — > e + e~ and (b) J/^/j — > ^ + /J-~ candidates. The 
dashed line in (a) shows the mass spectrum before the addition of bremsstrahlung photons. The 
arrows delimit the J/ip candidate region. 



from photon pairs detected within the barrel region of the Csl calorimeter, | cos# 7 |< 0.71, 
where 9 1 is the polar angle of the candidate photon with respect to the e + e~ beam axis. 
The photons must have a minimum energy of 30 MeV and their normalized invariant mass is 
required to be within 2.5a of the 7r° mass. This diphoton invariant mass is then kinematically 
constrained to the tt° mass. Charged and neutral pions and kaons are used to reconstruct 
the four K* decay modes. Candidate K* mesons are required to have a invariant mass 
within 75 MeV of the world-average K* mass [|J. 

We fully reconstruct 5-meson candidates by employing the kinematics of a BB pair 
produced almost at rest. We use the energy difference AE = E(J/ip) + E(K^) — -Ebcam as 
well as the beam-constrained mass M(B) = \J E^ cam — p 2 (B) as selection observables. The 
resolution in AE is 15 MeV for J/ipK* with a 7r° candidate in the final state and 9—11 MeV 
for the other modes. We find the resolution in M(B) to be ~2.5 MeV, which is dominated 
by the beam energy spread. We select signal candidates by requiring 5.2 GeV< M(B) < 
5.3 GeV and | AE |< 3cr ae- The beam-constrained mass distributions for events within the 
AE signal region are shown in Fig. F2|. 

We extract the signal yield in each mode by performing a binned maximum-likelihood 
fit to the M(B) projection, where the signal is given by a single Gaussian distribution with 
fixed mean of 5.28 GeV and fixed width of 2.5 MeV. The background is fit to a first-order 
polynomial joined with an elliptic function to fit the threshold nature of the beam-constrained 
mass distribution. The M(B) distributions in the AE sideband regions exhibit a slope 
consistent with zero. These sideband regions are at least 4 <tae and less than one pion mass 
away from the AE signal region. We fix the slope of the background shape to zero and 
allow the level of the combinatoric background to be determined from the fit to the M(B) 
projection of the AE signal region. 

We must account for the individual final states being reconstructed in a different channel 
(cross-feed), since for such candidates both the total energy and the beam-constrained mass 
lie near the signal region. We evaluate the reconstruction efficiency, as well as the amount 
of cross-feed from a given channel i to another channel j, using a sample of simulated B — > 
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FIG. 2. Beam-constrained mass projections (histograms) for candidates in the AE signal region 
are shown for the entire data sample summed over both J ftp — > l + l~ modes. The fits to the data are 
shown with the solid curves while the background fits are given with the dashed curves. Shown are 
distributions for (a) B+ -► J/^jK + , (b) B° -» J/?pK%, (c) B+ -► J/^jK*+, and (d) £° -> J/^if* 
candidates. 
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J I ipK^ events to generate a 6 x 6 efficiency matrix for the J/ip — > e + e and J/ip — > 
cases, as well as for CLEO II and CLEO II. V, separately. The CLEO detector simulation is 



based upon GEANT [JL2) . Simulated events are processed in a manner similar to that for the 
data. There is negligible cross-feed between the J/ipK and the J/ipK* modes. The cross- 
feed into J/ipK* modes with a charged-pion K* daughter is near 5%, whereas cross-feed into 
J/ipK* modes with a neutral-pion K* daughter ranges between 8 — 30% of the raw yield. 
Efficiencies and cross-feed-corrected yields are listed in Table |. As a cross check, we also 
quote the branching fractions computed for the analyzed B — > J/ipK^*' decays. 





efficiency [%] 


cross-feed corrected yield 


B 




J /if) — > e + e 


J/ip -> fj+fj, 


J/ip — ► e + e 


J/ip -> fi+fj, 


[xlO 






CLEO II 






J/^K + 


46.0 ±0.7 


56.2 ±0.7 


87.6 ± 9.4 ± 1.3 


121.9 ± 11.3 ± 1.5 


1.02 ± 0.07 


Jj^K% 


43.3 ±0.6 


51.9 ±0.7 


32.9 ±5.7 ±0.5 


24.0 ± 5.1 ±0.3 


0.83 ±0.12 


J/ipK*+(K^ir + ) 


28.0 ±0.6 


34.1 ±0.7 


13.1 ±3.8 ±0.3 


15.9 ± 4.6 ±0.3 


1.02 ±0.21 


J/^K* + (K+ir ) 


16.5 ±0.4 


20.7 ±0.4 


17.1 ±4.7 ±0.7 


23.6 ± 6.2 ±0.5 


1.63 ±0.31 


J/ipK*°(K+v -) 


33.5 ±0.5 


40.6 ±0.5 


53.1 ±7.5 ±0.8 


56.8 ± 8.0 ±0.7 


1.11 ±0.11 


J/i/jK*°(K° s ir ) 


16.4 ±0.4 


18.7 ±0.4 


3.4 ±2.3 ±0.2 


5.1 ± 3.0 ±0.2 


1.02 ± 0.46 






CLEO II. V 






J/ipK+ 


43.6 ±0.7 


58.8 ±0.8 


172.5 ± 13.1 ±2.8 


210.2 ±15.0 ±2.9 


0.98 ± 0.05 


J/%bK% 


42.5 ±0.5 


57.9 ±0.6 


42.5 ± 6.5 ±0.6 


78.9 ± 9.1 ±0.9 


0.90 ± 0.08 


J/ipK*+{K° s iT + ) 


28.6 ±0.6 


36.4 ±0.7 


17.0 ± 4.6 ±0.5 


47.6 ± 7.3 ±0.9 


1.00 ±0.14 


J/lpK*+{K+TT°) 


14.8 ±0.4 


21.6 ±0.4 


35.9 ± 6.8 ±1.0 


42.0 ± 8.2 ±1.0 


1.69 ±0.23 


J/ipK*°(K+ir ) 


31.1 ±0.7 


41.1 ±0.8 


92.5 ± 9.8 ±2.1 


105.8 ±11.1 ±2.1 


1.08 ± 0.08 


J/^K*°(K s tt°) 


15.5 ±0.5 


18.2 ±0.5 


11.9 ± 4.0 ±0.5 


8.6 ± 4.0 ±0.4 


1.37 ±0.38 



TABLE I. Summary of reconstruction efficiencies (daughter branching fractions not included), 
cross- feed corrected signal yields (first error is statistical, second error is systematic), and branch- 
ing fractions B computed from these yields and efficiencies (errors are statistical only), assuming 
equal production of B + B~ and B°B pairs. The computed branching fractions agree with the 
world-average values Q. Results for the data accumulated with the CLEO II and CLEO II. V 
configurations are given separately. 

We extract our result from the cross-feed and reconstruction-efficiency corrected yields 
using the world-average values for the respective daughter decay branching fractions []]]. We 
obtain four independent measurements of R listed in Table 

We evaluate the uncertainties in the reconstruction efficiency due to track finding, track 
fitting, charged hadron identification, Kg finding, and tt° finding. Since we use the ra- 
tio of two decay rates that each involve J/ip — > l + l~ candidates, uncertainties in lepton 
identification are negligible. We estimate the full systematic bias due to daughter recon- 
struction efficiency uncertainties by taking into account correlations between the different 
final states in the numerator and denominator, resulting in some cancellation. Propagating 
these uncertainties through the weighted average of the results in Table H, we arrive at a 
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Configuration 








CLEO II 


CLEO II. V 


Signal Mode 






B -» J/ipK 


1.229 ±0.191 


1.088 ±0.116 


B -» J/V'K* 


1.098 ±0.190 


1.095 ±0.137 



TABLE II. Results for R = 4t- x for the different CLEO configurations and the JlibK 

TOO T B o o I Y 

and J/tpK* modes. The uncertainties are statistical. 

systematic uncertainty on R due to the understanding of reconstruction efficiencies of ^15%. 
The polarization of the decay B — > J/ipK* is modeled in our simulation with a longitudinal 
polarization fraction of Yl/Y = 0.52 in accordance with Ref. 0. We estimate the impact 
of the value used for Yl/Y on the B — > J/ipK* efficiencies by generating signal events with 
the nominal polarization varied by ±lcr = ±0.08. The central value for R changes by less 
than 0.8% due to this variation. We vary the B candidate signal width by ±0.2 MeV and 
estimate the systematic uncertainty from this source to be less than 0.5%. We extract the 
signal yield using a background function that allows for a slope in the non-signal region of the 
beam-constrained mass distribution and assign a systematic bias of 3% on the central value 
for R from our assumption of flat background. We attribute a 1.1% uncertainty to limited 
statistics of the simulated event samples used to extract the efficiency matrices. Adding all 
contributions in quadrature we obtain a total systematic uncertainty of ^37% on R. 

We weight the results of Table [Tj with their statistical uncertainty and, combined with 
the estimated systematic uncertainty, we extract 

R = I±- x Zs+ = 1.11 ± 0.07 ± 0.04, 

700 t s o 

where the first uncertainty is statistical and the second is systematic. 
Using the world-average lifetime ratio of charged and neutral B mesons, 1.066±0.024 [13 
we obtain a measurement of the production ratio 



= r(T(4S)^B+B ) = LQ4 ± Q Q7 ± Q 
Too r(T(4S)->B°B ) 

and, assuming / + _ ± / o = 1? we also extract /oo = 0.49 ± 0.02 ± 0.01 and /+_ = 
0.51 ±0.02 ±0.01. 

We have measured the ratio of charged to neutral production of B mesons at the T(AS) 
resonance to be 1.04 ± 0.07 ± 0.04, which is consistent with unity within an error of 
8%. This is the most precise measurement of /+-//oo- Our result is also consistent with 
theoretical predictions of greater charged than neutral 5-meson production in T(4S') decays 
near threshold. We emphasize that the ratio / + _//oo and its uncertainty must be taken into 
account when performing measurements that compare charged and neutral B decays at the 
T(4S') resonance. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences and 
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